Fluctuating Nematic Elastomer Membranes: a New Universality Class 
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We study the flat phase of nematic elastomer membranes with rotational symmetry spontaneously 
broken by in-plane nematic order. Such state is characterized by a vanishing elastic modulus for sim- 
ple shear and soft transverse phonons. At harmonic level, in-plane orientational (nematic) order is 
stable to thermal fluctuations, that lead to short-range in-plane translational (phonon) correlations. 
To treat thermal fluctuations and relevant elastic nonlinearities, we introduce two generalizations 
of two-dimensional membranes in a three dimensional space to arbitrary D- dimensional membranes 
embedded in a d- dimensional space, and analyze their anomalous elasticities in an expansion about 
D = 4. We find a new stable fixed point, that controls long-scale properties of nematic elastomer 
membranes. It is characterized by singular in-plane elastic moduli that vanish as a power-law 
r/A = 4 — D of a relevant inverse length scale (e.g., wavevector) and a finite bending rigidity. Our 
predictions are asymptotically exact near 4 dimensions. 
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I. INTRODUCTION 

The ubiquity and importance of membrane realizations 
in nature, as cellular walls, and in laboratories, as self- 
assembled bilayers of lipid amphiphils, has fueled consid- 
erable scientific activities U. These nearly tensionlessQ 
lipid sheets are highly flexible, with elastic moduli of- 
ten comparable to thermal energies. Consequently, on 
long length scales, their conformational properties are 
strongly affected by thermal fluctuations. This, together 
with early successes in understanding a variety of puz- 
zling phenomena (such as for example red blood cells' 
flicker,^, biconcave shape of eurythrocites^l, and period 
of lyotropics) in terms of continuum models of fluctuating 
elastic sheets, has attracted the attention of the physics 
community. Consequently, significant progress has been 
made in understanding the statistical mechanics of fluc- 
tuating membranes 

It is by now well-appreciated that the nature of a mem- 
brane's in-plane order, with three (heretofore studied) 
universality classes, the isotropic, hexatic and solid (i.e., 
tethered or polymerized), crucially affects its conforma- 
tional properties. The most striking effect of in-plane 
orders is the stabilization in solid membranes of a "flat" 
phase 5] , with long-range orientational order in the local 
membrane normals 6] , that is favored at low-temperature 
over the entropically preferred high-temperature crum- 
pled state. Therefore, in marked contrast to liquid 
membranes and one-dimensional polymer analogs, which 
are always crumpled (beyond a persistence length) jTJ, 
tethered membranes, despite of being two-dimensional [g 
are predicted 5] to undergo a thermodynamically sharp 
crumpled-to-flat phase transition 9, 10]. The ordering is 
made possible by a subtle interplay of thermal fluctua- 
tions with nonlinear membrane elasticity, which, at long 
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scales infinitely enhances a membrane's bending rigid- 
ity, thereby stabilizing the orientational order against 
these very fluctuations. This novel "order from disor- 
der" phenomenon, and the universal "anomalous elastic- 
ity" , namely, length-scale dependent elastic moduli, non- 
Hookean stress-strain relation, and a universal negative 
Poisson ratio 0, 0, 0, are now known to be quite 
commonly exhibited by many other "soft" systems sub- 
jected to fluctuations. 

In this paper we introduce and explore a new univer- 
sality class of solid membranes, that spontaneously de- 
velop an in-plane orientational nematic order. Our mo- 
tivation is two-fold. Firstly, our interest is driven by 
experimental pro gres s in the synthesis of nematic liquid 
crystal elastomers TS], statistically isotropic and homoge- 
neous gels of crosslinked polymers (rubber), with main- 
or side-chain mesogens, that can spontaneously develop 
nematic orientational order. Even in the absence of fluc- 
tuations, they were predicted^ and later observed to 
display an array of fascinating phenomena|l5j . the most 
striking of which is the vanishing of stress for a range 
of strain, applied transversely to the nematic direction. 
This striking softness is generic, stemming from the spon- 
taneous orientational symmetry breaking by the nematic 
state0, 01 that ensures the presence of a zero-energy 
Goldstone mode, corresponding to the observed'2l5 soft 
distortion and strain-induced director reorientation. The 
hidden rotational symmetry also guarantees a vanish- 
ing of one of the five elastic constants 11^ that usu- 
ally characterize harmonic deformations of a three- 
dimensional uniaxial solid 22] . Thermal fluctuations lead 
to Grinstein-Pelcovits-like^lTj renormalization of elastic 
constants 01 in bulk systems with dimensions below 3 
in pure systems 0, 0| and below 5 when effects of the 
random network heterogeneity is taken into account 23|. 
It is, therefore, likely, and indeed we find, that the elas- 
tic properties of a two-dimensional fluctuating tcnsionless 
sheet of such nematic elastomer differ qualitatively from 
those of previously studied crystalline membranes Q. 
Our aim here is to explore the effects of thermal fluctua- 



2 



tions on this new universality class of solid membranes. 

Our other motivation for exploring the physics of ne- 
matic elastomer membranes comes from earlier discov- 
ery by Toner and one of us|23j|. that any amount of 
any kind of in-plane anisotropy, a seemingly innocuous 
generalization, significantly enriches the phase diagram 
of polymerized membranes. Most dramatically, it was 
predicted (2il| that an entire new phase of membranes, 
called the "tubule" phase, always intervenes between the 
high-temperature crumpled and low-temperature "flat" 
phases. The defining property of the tubule phase is that 
it is crumpled in one of the two membranes directions but 
"flat" (i.e., extended) in the other, with its thickness, 
roughness, and anomalous elasticity displaying univer- 
sal behaviors controlled by a nontrivial, infrared stable 
fixed point The tubule phase has since been observed 
in Monte Carlo simulations {24| of non-self-avoiding (i.e., 
phantom) anisotropic membranes with properties (the 
thickness, roughness, and Poisson ratio) in close qual- 
itative and quantitative accord with the predictions of 
Radzihovsky and Toner '23'|. 

Although it was an explicit rotational symmetry break- 
ing (ariisotropy) that was considered by the authors of 
Ref. !23j we expect a similar, but qualitatively distinct 
tubule phase to be also displayed by the spontaneously 
anisotropic elastomer membranes considered in this pa- 
per. However, we leave the subject of the global phase 
diagram, the elastomer tubule phase, and the phase tran- 
sitions between it and the flat and crumpled phases for 
a future publication psf. Here we instead focus on the 
formulation of models of nematically-ordered elastomer 
membranes and the study of long-length-scale properties 
of its flat phase. 

Past extensive investigations, both in the context of 
defected crystalline membranes '26^, '27] and liquid crys- 
tals confined in rigid gels (e.g., aerogels and aerosils) 28J, 
demonstrated that arbitrarily weak heterogeneity qual- 
itatively modifies the long-scale nature of liquid-crystal 
orders, membranes morphologies and elasticities. We, 
therefore, expect that local heterogeneity in an elastomer 
network will also become qualitatively important on suffi- 
ciently long scales and will likely dominate over the ther- 
mal fluctuation effects that are the subject of this paper. 
However, here we consider the idealized limit of elastomer 
membranes in which the effective quenched disorder is 
sufficiently weak that its effects are unimportant on scales 
shorter than some very long disorder-determined length 
scale. We leave the study of these heterogeneity effects 
for a future investigation|29| . 

This paper is organized as follows. In Sec. m we de- 
velop two models for /^-dimensional nematically-ordered 
flat phase of elastomer membranes fluctuating in d- 
dimensions with d > D, one with uniaxial nematic or- 
der, the other with D-axial nematic order. They are re- 
ferred as uniaxial model and D-axial model respectively 
in this paper. For the physical case of interest D — 2, 
they reduce to the same model, which we study in more 
detail in Sec. IIVI In Sec. IIIII we investigate harmonic 



fluctuations of both models and show, in particular, that 
in-plane phonon fluctuations in D > 3 have parts that 
are equivalent to those of smectic liquid crystal, a colum- 
nar liquid crystal, and a crystalline solid. In Sec. IIVI we 
investigate mean-fleld theory and harmonic fluctuations 
in physically realizable two-dimensional membranes em- 
bedded in a three-dimensional space. We flnd that har- 
monic fluctuations in these membranes do not destroy 
long-range in-plane nematic order in spite of enhanced 
fluctuations relative to an isotropic system arising from 
the soft Goldstone mode. In Sec. we consider the 
effects of anharmonicities and develop an e = (4 — D)- 
expansion about the upper critical dimension, Dc — 4, for 
two model systems that have well-deflned analytic contin- 
uations to D = 2. In both models, we find similar flow 
equation for coupling constants and an infrared-stable 
fixed point in which the bending modulus k is unren- 
mormalized and in-plane elastic moduli A,., vanish with 
wavenumber q as q'^. 



II. MODELS 
A. Notation: Reference and Target Spaces 

Although physical membranes arc two-dimensional 
manifolds in a three-dimensional space, it is often useful 
to consider generalizations to _D-dimensional manifolds 
embedded in a d-dimensional space with d > D. Our 
interest is in solid membranes, which, unlike fluid mem- 
branes, have a non-vanishing shear modulus at least in 
their isotropic state. 

To describe the geometry and fluctuations of these 
membranes, we need to introduce a certain amount of 
notation. First, we deflne the reference space Sr. This is 
the £)-dimensional space occupied by the membrane in its 
quiescent flat reference state. We denote D-dimensional 
vectors in this space in bold and their components with 
Roman subscripts i,j = 1,...,D . In particular, we de- 
note intrinsic coordinates on the membrane by the vector 

r ^ {xi, ...,xd) (2.1) 

with components Xi. The positions, r index mass points 
in a Lagrangian description of the membrane. The mem- 
brane fluctuates in a c?-dimensional target space St- 
We denote vectors in the embedding space with arrows 
and components of these vectors with Greek subscripts 
fijV — 1, ...,d. In particular, we describe the position in 
St of the mass point labeled by r with the embedding 
vector 

i?(r) = (i?i(r ),..., i?,(r)) (2.2) 

with components -R^(r). 

The reference space Sr is identified as a subspace of 
St- The use of orthonormal basis vectors will simplify 
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some of our discussion, and we introduce d vectors in 



St with components e^^ — 6^^^, satisfying 



(2.3) 



Any vector in St can be decomposed into its components 
along these vectors. For example R = Rfj,e^, where the 
summation convention on repeated indices is understood. 
We choose the set {e^} so that the subset of vectors e^, 
i = \, ...,D lie in Sr, which we represent as a subspace 
of S't. Thus 



1 iii = fi<D 
if = £» + 1, 



(2.4) 
(2.5) 



u, with components Uy = gij — glf , measures the local 
change in separations of nearby points in distorted states 
characterized by ^(r) relative to those of the undistorted, 
flat state characterized by i?rof: 



dR^ — dR^„f — 2uiidxidxj . 



(2.12) 



Using Eq. H2.7|l . the strain can be expressed in terms of 
the phonon and height fields as 




OR dR 



diU ■ djU + dih ■ djh)(2.13) 



We choose length scales such that the position in St of 
the point r in the reference membrane is 



Rvcf — XiCi 



(2.6) 



describing a perfectly flat membrane configuration. Po- 
sitions R{r) of the distorted membrane can be described 
by a Z?-dimensional in-plane phonon field u(r) and a 
d — D = o?c-dimensional out-plane undulation (height) 
field h{r): 

D 

R{r) = ^{x^ + u,(r))e, + ^ hk{r)ek+D- (2.7) 
1=1 fc=i 

Local distortions of the membrane can be expressed in 
terms of the d x ZJ-diniensional deformation tensor A j^^l 
with components 



A„ 



dR^ 
dxi 



In the flat reference state, 

of _ dRici^fj, 



Ay: 



dxi 



(2.8) 



(2.9) 



Under independent rotations in the target and reference 
spaces, Api transforms as 



A 



(2.10) 



where 0T/ii/ is a d x d rotation matrix in St and 0_R,y a 
D X D rotation matrix in Sr. 

Since i? is a Euclidean vector, the metric tensor, 
which measures distances between neighboring points via 



9ij 



dRjr) dRjr) 

dxi 



dx-i 



In the reference flat phase, gij = gl^ — A^'AJ^'^j = Sij. 
We will denote tensors in Sji with a double underscore. 
Thus the metric tensor is g. The Lagrangian strain tensor 



B. The Model Hamiltonian 

As discussed in the Introduction, internally isotropic 
solid membranes have been extensively studiedjlj. In 
contrast to liquid membranes and one-dimensional poly- 
mer analogs, they admit a flat phase, characterized by 
an infinite orientational (in embedding space) persistence 
length that is stable to thermal fiuctuations. The long- 
wavelength free-energy density for distortions of such 
spontaneously flat phase is given by 



where 



and 



= -^bend + -^stretch, (2-14) 



^bc„d = :r(v2;i)2 (2.15) 



stretch = ^{T^Tgf + fiiTm^) (2.16a) 
B 



{TtuY + Ai(TrM^), (2.16b) 



where k is membrane's bending rigidity A and fi are 
the Lame coefficients characterizing in-plane elasticity, 
B = X + 2fi/D is the bulk elastic modulus, and 



— (Tr^)5,_ 



(2.17) 



is the symmetric-traceless part of the strain tensor with 
Tru — 0. Earlier studies of non-liquid crystalline poly- 
merized membranes 1^ 0, Il2| have demonstrated that 
thermal fluctuations produce wild undulations about the 
flat state Rref that make elastic nonlinearities in the 
height field h (but not in the in-plane phonon fields 
(2 11) u) in Uij Eq. H2.13|l important on long length scales 
and lead to universal length-scale-dependent elastic mod- 
uh: K(q) g"''", A(q) ~ q^i^ , and ^(q) ~ g''". One 
of the most important consequences of this is that the 
thermally-driven upward renormalization of the bend- 
ing rigidity k stabilizes the low-temperature flat phase of 
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two dimensional H membranes against these very same 
fluctuations 0, 0, . 

The free energy density of Eq. H2.14|l provides a cor- 
rect description of elastic and height fluctuations of a 
membrane so long as the equilibrium phase is truly an 
isotropic flat phase. If, however, the shear modulus fi 
becomes sufficiently small, the membrane becomes un- 
stable to spontaneous in-plane distortions that break ro- 
tational symmetry of 5'_r0, 0, ll^. In the presence 
of such instability the distortion is stabilized by a non- 
linear strain energy ^jvl and an in-plane curvature en- 
ergy .7-"curv, neglected in JF/, Eas l2.14l2.1t)bl The in-plane 
curvature energy, which stabilizes the system against spa- 
tially nonuniform distortions is simply 



2„\2 



^curv - -X(V"U 



(2.18) 



The important low-order contributions to non-linear 
strain energy are 

Tnl = ^TtuTtu^ - CTr(u)3 + "fiTiu^. (2.19) 



specified by a unit vector n° with components n?, is ar- 
bitrary in the plane Sn, and we could take it without loss 
of generality to be along one of the basis vectors, say ii. 
The vector n° also exists in the target space St- It has 
components nj^ that are equal to for fj, — i — 1, _D 
and zero for fi > D. Thus, in the uniaxial case, 

Ay- = (Ao|| - Ao±)noinoj + Ao±S^J. (2.22) 

Under rotations in St and Sj^, this becomes 



(2.23) 



where we used S^i — e,. 



and where 



= 07 



= &T^,^e^^, ej)^ = Qfj^^j, and nf = QfjUoj. 
The opposite limit of the uniaxial anisotropic state is 
the full Z3-axial state in which there is unequal stretching 
in all D directions in the plane. In this case, 



D 



(2.24) 



Terms of order or higher and terms proportional to 
(Trw)^ and (Tru)^Tru^ could also be added but are qual- 
itatively unessential for our current discussion. 

C. Spontaneously Anisotropic Phases 



When fi becomes sufficiently small in the presence of 
Tnl, a membrane will undergo a transitionp^ IT^ l25l | 
to a new anisotropic, spontaneously stretched equilib- 
rium state with a non- vanishing equilibrium strain u^, 
an anisotropic deformation tensor Aj^j that differs from 
A)j'i^ = 6^i, and position vectors i?o,Ai = A^i^i- 



The 



most general form of A° ^ can be obtained via target- and 
reference-space rotations Eq. H2.10|) of the deformation 
tensor 



A' 



ii fi = i ~ 1, D 
ii fi = D + l,...,d 



(2.20) 



resulting from a distortion without rotation in the Sr 
subspace of St- Here A°^ are the components of the a 
D X D matrix A° restricted to the plane of the origi- 
nal reference space Sb,- The new equilibrium A" yields 
a new metric tensor g^j = AjJ^A^^ and a non- vanishing 
equilibrium strain relative to the original isotropic flat 
membrane: 



with all Aofe different. 

As discussed in detail in Ref . , the spontaneous bro- 
ken symmetry described by the above deformation ten- 
sors leads to vanishing of certain shear moduli of the dis- 
torted solid. Distortions relative to the new anisotropic 
state are measured by the displacement and height vec- 
tors, u'(r') and /i(r') defined via 



r' + u'{r') + h{r')=R{r), 



(2.25) 



where by definition, r' — Rq, u'(r') is the the D- 
dimensional in-plane displacement vector, and h{r) is the 
dc-dimensional {dc = D — d) height variable orthogonal 
to u'. The strain tensor relative to the new state, 



, _ 1 dR^dR^ . 
" 2 dx' dx', 



(2.26) 



expressed in terms of the strain tensor u relative to the 
original state is 



(2.27) 



The strain Hamiltonian density to harmonic order for the 
uniaxial case 



-'^strain ~ r, / , "^ij /^-L 



(2.28) 



(2.21) 



The simplest anisotropic state that can form is the 
uniaxial one in which the original isotropic membrane is 
stretched (compressed) along a single direction in Sn and 
compressed (stretched) along the others. Because the 
reference state is isotropic, the direction of stretching. 



where the Einstein convention is not used in the first term 
on the right hand side (but is in the second). 



3y - riojnoj. 



(2.29) 



(2.30) 
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and 



X-f = Ai5,i<5,i + A3(l~fe)(l"<5ji) 

+A2[(5a(l - 6ji) + (5,i(l - da)], (2.31) 

where the values of Ai, A2 and A3 depend on the po- 
tentials of the original Hamiltonian and Aq. This form 
makes it clear that Ay is actually a subset of a fourth rank 
tensor and not a true second-rank tensor. The bending 
and curvature energies become anisotropic with 



direction, perpendicular to the direction of order, and 
^ij" = ^iy^jy'^'yy Thus the second term in Eq. (|2.28|) can 
be absorbed into the first term, which depends only on 
diagonal components of Uij. Thus, in 2D, the free energy 
density has the same form as Eq. H2.35|l . and is most nat- 
urally analytically continued from _D-axial generalization 
of a nematic elastomer membrane. 



III. HARMONIC THEORY 
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where the anisotropic bending and curvature moduli can 
be expressed, respectively, as Ky = k^™' and Kij = 
KA^'f where 

''J 

Aj;- = A^||<5,i^,i + A4^(1 - 6a){l - 5,,) (2.34) 
- + <5ji(l - 5a)]- 

In the full Z3-axial case, there is no residual plane that 
can support shears and the harmonic free-energy density 
depends only on the diagonal parts of the strain: 



-r-D-axial _ }_ \ i i 
•' strain n "ii "jj 



(2.32) The two (uniaxial and D-axial) nematic membrane 
Hamiltonians introduced in the preceding section has 

(2.33) complicated anisotropics, associated with spontaneous 
in-plane nematic order, that lead to membrane's highly 
anisotropic conformational correlations. In this section, 
we will investigate fluctuations in the harmonic approx- 
imation to these models. These harmonic results will be 
necessary in our subsequent discussions in Sec. [V] of the 
important effects of nonlinearities in the presence of fluc- 
tuations. To simplify the notation, in this section we will 
drop primes appearing in the strain and phonon fields 
and u' measured relative to the nematic state and denote 
them by Uij and u, not to be confused with physically 
distinct unprimed quantities of the previous section. 

The harmonic Hamiltonians for the uniaxial and D- 
axial models have a longitudinal-strain and height elastic 
(2.35) parts describing the cost of simple shear modes: 



with Ajj a symmetric matrix with all D{D + l)/2 inde- 
pendent entries different. Again, the Einstein convention 
is suspended in this equation. The bend and curvature 
energies have the same form as Eqs. (|2.33|l but A"^"' is 
replaced by 



1 f d^q 



(3.1) 



i=D+l 



4 D — axial a2 a2 



(2.36) where 



We note that in the most general case that is compatible 
with D-axial symmetry, all D{D + l)/2 components of 
^D-axiai g^^g independent of each other. 

The total Hamiltonian density for a nematic elas- 
tomeric membrane is 



K{q) = Y^K^.qfq 

ij 



(3.2) 
(3.3) 



(2.37) 



We have expressed .Fstrain only to harmonic order in the 
nonlinear strain Uij even though nonlinear terms in the 
original isotropic energy are essential to stabilize the sys- 
tem after nematic order develops H^. As we shall see in 
Sec. O near four dimensions anharmonic terms in the 
nonlinear strains associated with in-plane phonons are 
subdominant (irrelevant in RG sense) to nonlinearities 
in the height undulations h and we will ignore them in 
what follows. This is in strong contrast to bulk nematic 
elastomers, where such phonon nonlinearities cannot be 
ignored and must be treated non-perturbativelv |l8l IT9| . 

In two dimensions, which we consider in more detail 
in the Sec. IIVI there is only one direction, say the y 



with qi — qil q. In uniaxial systems, there is in addition a 
shear part in the (13 — 1) dimensional plane perpendicular 
to the unique direction Hq. 



(3.4) 



where is the part of the vector u in the plane perpen- 
dicular to Hq. u-'- has D — 1 components. Thus in Z3 = 2, 
has only one component, say Uy, ujj = Uyy and two 
models are equivalent. 

In both models, at harmonic level, height and phonon 
variables decouple, and the height correlation function is 
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simply 



Ghihj — Sl 



(3.5) 



di- 



where Sfj is the projection operator onto the dc 
mensional subspace perpendicular to the reference mem- 
brane. 



A. Uniaxial Case 

In the uniaxial case, u can be decomposed into compo- 
nents along the uniaxial direction no and along directions 
parallel and transverse to the wave- vector qj^ in the plane 
perpendicular to np: 



UnUOi + ULq±i + Uti, 



(3.6) 



where no • Ut = and qj^ ■ Ut = 0. We will represent this 
decomposition as Ui = (M„,u;,Ut). The u-correlation 
function can then be decomposed as 

Gij = GnnnoiriQj + GLLq±iq±j + GnL{nQiq±j + nQjq±i) 
+Gt{5l-q^,q^,). (3.7) 

The transverse part of decouples from the others and 
is easily calculated 



Gt 



1 



(3.8) 



The Hamiltonian for the Un and ul components can be 
written as 



where u 



G- 



n 



and 



Aig2 + K{q)q^ 
A2'7||g_L 



A2g||g_L 
\;ql + K{q)q^ 



(3.9) 



(3.10) 



where A3 = A3 + /ij^ with Aj defined in Eq. (|2.31|l . In the 
long- wavelength limit, terms of order gj^ and q^\q\ can be 
neglected relative to gj^ , and the nn entry in G~^ is effec- 
tively Aigj^H-Xj^g^, which is the inverse of the propagator 
for fluctuations in a smectic liquid crystal. Similarly, in 
the long- wavelength limit, the LL part of G~^ is effec- 
tively \^q\ + which is the inverse propagator for 
one component of displacement fluctuations in a colum- 
nar liquid crystal. Thus, when A2 = 0, G^^ decomposes 
into two independent parts, one of which is identical to 
the inverse propagator of a smectic liquid crystal and the 
other of which is identical to one component of the prop- 
agator of a columnar liquid crystal. We will find below 
that this property survives the turning on of A2 though 
with different coefficients. 



Equation H3.10|l is easily inverted to yield 
1 / A3a^ +K(q)q'' -X 



G = 



A 



hqi+ K{q)q^ -A2q\\qj_ 
-X2q\\q± Xiqf^+ K{q)q'^ 



(3.11) 



where 

A = dct G-' = A^gig2 + (Aig2 + X3ql)K{q)q^ + ■■■ 

(3.12) 

with A;^ — A1A3 — A2. It is can be shown (See Appendix 
A for discussion in two dimension) that Gnn is dominated 
at small q by the region with ^ g^ so that 



G^-j/rt, 



Xml + K{q)q^ 



where Ai = A^/Ag = Ai - A2/A3, Ky_ = L(Ao_l 
larly, Gll is dominated by q±_ « g^ at small q and 



(3.13) 
Simi- 



Gll — 



Aig^ + K{q)q^ 

A 



(3.14) 



where X'^ = A-JXi, K^^ ^ K(Aoii)\ 

Thus we find, that, within the harmonic approxima- 
tion, the translational lower critical dimension is de- 
termined by the dominant smectic-like fluctuations of 
Un modes and is equal by 13 = 3. The subdominant 
columnar-like modes become important only below 
D — 5/2, while the transverse fluctuations Ut only grow 
with length scale for D < 2. 



B. D-axial systems 

In Z3-axial systems, the strain energy is given by Eq. 
(|2.35() with D{D + l)/2 independent components for Xij. 
The harmonic Hamiltonian is then given by Eq. 13. 1|) 
with no transverse contributions. The height fluctuations 
are the same as in the uniaxial case with the appropriate 
expression for K{q). The inverse phonon propagator is 



^Xiql + K{q)q^ 



MoqiqD 



Aingito 



Xoqh + K{q)q^ 



(3.15) 



G can be obtained by inverting G ^ using minors: Gij = 
(-ly+JA^yA where A = detG^^ and A*^' is the deter- 
minant of the minor of G~^ obtained by deleting row i 
and column j. The results for the long- wavelength diago- 
nal and off-diagonal components of Gij are, respectively. 



Gi 

G., 



Xyf + Kiq)q^^ 

{~iy+^A=^q.q, 



^xq 



2„2 



q.+K{q)q^A^^q: 



(3.16a) 
(3.16b) 
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where, as before, A a is the determinant of A, is the 
determinant of the ij-minor of A, and AJj = det Aa/A^-*. 
Thus, the diagonal components Gu have the structure 
of a smectic propagator, with all terms involving qi in 
K{q)q^ subdominant in the long wavelength limit. 



IV. TWO-DIMENSIONAL MEMBRANES 

As discussed in the Introduction, real membranes are 
two-dimensional manifolds, with reference-space coordi- 
nates r — {xi,X2)^ fluctuating in a three-dimensional 
space with coordinates R — {Ri, R2, Rs). Because of 
their low dimensionality, these membranes have a number 
of special properties, some of which we will investigate in 
this section. In particular, we will study the transition 
from an isotropic to a nematically-ordered membrane, 
which, because of the absence of a third-order invariant 
for two-dimensional symmetric-traceless tensors is of sec- 
ond order in mean-field theory. Our treatment will re- 
view in a two-dimensional context the mechanism for the 
emergence of soft elasticity in the nematic phase. We 
will also consider harmonic fluctuations about the flat 
nematic state in some detail and show that they lead to 
short-ranged positional correlations, but allow for a sta- 
ble long-range orientational (nematic) in-plane order of 
the membrane. 



value Tru^ at equilibrium is zero so long as /i > 0. In 

this case, Tru^ is also zero. When /i < 0, Tru^ and Tru^ 
become nonzero with values 



B\f,\ 



Tru? 



(3 2 
Trio = -^Tr^o' 



2B7 - /32 ' 



(4.2a) 



(4.2b) 



2S7 - /32 ■ 

By choosing axes so that the anisotropic stretching is 
along the x-axis, we can express and in terms of 
Aoa; and Aoy as 

(4.3b) 



A2^ - 1 



i^Ox ^Oy) ( — 1 



From this and 



u = — (Tru )I + u 
=0 2 =0'= =0 



1 



(4.4a) 



we obtain 



The Two-Dimensional Hamiltonian 



Ao, = ( 1 + Tr|^ + V2Tr£2 j , (4.5) 



Our starting point is the elastic energy of an isotropic 
two-dimensional membrane augmented by nonlinear elas- 
tic terms necessary to stabilize the nematic state that 
develops when the shear modulus becomes negative. A 
simplifying feature special to two dimensions is that the 
strain tensor u is a 2 x 2 matrix with no independent 
cubic invariant. The corresponding Hamiltonian is 



|(V2/I)2 + f(V^u)2 + |(Tr^)2 



/3TruTru^ + 7(TrM 



2 

(7 



2B' 



f(V^u)^ + |(Tr^- 



|Trr) 



(4.1a) 



,-,2 



m/2 



7 - PVim 



(4.1b) 



To simplify our discussion, we have left out the cubic and 
quartic terms, (Tr^)^, (Tr^)"* and (Tr^)^(Tr£^), as well 
as higher order nonlinearities in u. The inclusion of these 
terms, which are quantitatively smaller than the terms 
we keep in the nearly incompressible (large B) limit, will 
not qualitatively modify our results. 

The nonlinear strain tensor u'^j Eq. (|2.21|) associated 
with the spontaneous deformation in the nematic state 
can easily be found by minimizing the effective energy 
Eq. H4.1|l . Since Tru^ is greater than or equal to zero, its 



Aoy = fl + Tru 



2TtuI 



(4.6) 



where Tru^ and Tru^ are given in Eq. I|4.2|l . 

We can now massage the free energy Eq. H4.1|l , into 
a more convenient form, 



^ = -(V^/^)^ + -(V^u) 



+ y (Tr^ - Tr^o)^ + 7(Tr^'' " Tr£^)^ 

~,2 ^^„~,2^ {4:-7) 



+ (3{Ttu — Ttu^)(Ttu — TrWg), 

that makes it transparent that J- is minimized by u = 
and that permits a straightforward expansion about the 
ground state. 

Using Eq. H2.27|l . we can now easily express the free- 
energy density of Eq. 14.7|l in terms of the strain u' rel- 
ative to the new stretched equilibrium state. First we 
observe 



Ttu-Ttu = TrAV, 
= =0 =0 ' 

Tru^ _ Tru^ = 2Tr(u A^u') + Tr{A^u' f 
= =0 ^=0=0= ' ^=0= ' 



(4.8a) 



-(TrA^u' 
2^ =0= 



2„/^2 



(4.8b) 
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The final expression, valid to linear order in m', does not 
depend on u'^y - a property, whose origin is the sponta- 
neous broken rotational symmetry of the nematic phase, 
that is responsible for the vanishing of the membrane 
shear modulus. Using Eq. (|4.8|l in Eq. I|4.7|l . replacing di 
by dl = Aoijdj, and retaining only the dominant terms 
in 9,;U, we obtain 



1 



Md'yU, 



(4.9) 



where to streamline our notation we again dropped 
primes, by replacing u' with u and 9^' with di . The 
strains are the usual nonlinear strains relative to the 
new reference state, which to linear order are simply 



dxUx and Uy 



dyUy. 



The bare elastic coefficients 



^xy 

Kyy 

Kx 



^^Ox^^Oy'^7 
^ty^7 

AtyK, 



Ky = X, 



^tx(B 



-/3(a: 



2 

Ox 



Ky) 

Ky) 



>^xy = ALAL(S-7(A2^-Ai 



Oy) 



liKx 



Ky)'), 
Ky)% 



(4.10) 



B. Fluctuations and Correlations of the Harmonic 
Model 

We study fluctuations of nematically-ordered elas- 
tomer membranes within a harmonic approximation, for 
the physically realizable case of I? = 2 and d = 3. In this 
case, the displacement vector u has two components, and 
the height has a single component h. Within this ap- 
proximation all correlation functions are related to the 
harmonic two-point correlation functions |3l| 



Glir) 



G°-(r) = 



{h{r)hiO))o 



(2^)2 
Ui(r)uj(0))o 

{2ny 



e*''""G°,(q), 
e'''-""G°(q), 



(4.11a) 



(4.11b) 



expressed in terms of corresponding "propagators" Gf^{q) 
and G'-'j(q). As usual, the averages are computed using 

a Boltzmann weig ht Zo"^e-^«['''"l (for convenience using 



ksT as the energy unit), by integrating over phonon u 
and height undulation h fields, with Zq = J DhDu e^'^" 
the partition function and Tio the harmonic effective 
Hamiltonian Ho = J (PxJ-o[h,u] obtained from effective 
Hamiltonian, H — J (PxT[h, u] by neglecting all elastic 
nonlinearities appearing in Eq. (|2.13() . 

The height propagator G° (q) is given by Eq. (|3.5|l and 
the phonon propagator G°j(q) by Eq. (|3.16() specialized 
to two dimensions: 



GL(q) 

Gyy{(l) 

G"y(q) 



1 



Klx + Kylt' 



X'yQl + K^qt ' 



Axqlql + K{q)q\Xyql + X^q^) 



(4.12a) 
(4.12b) 
(4.12c) 



k's, K's, and A's are determined by the parameters of the where 
original Hamiltonian and A : 



X' — X 

'^x,y ~ ^x,y 



\2 

^xy 

At- A,, 



(4.13) 



>0, Xy> 0, and X^Xy-Xly > 0. A^ 
Xly. The domi- 



Stability requires A^: 
and X'y both go to zero linearly in Aj,Aj^ 
nant parts of the Ux and Uy correlation functions are iden- 
tical in form to the displacement correlation functions of 
a two-dimensional smectic, which have been extensively 
studiedjSO']. Both (u^) and (Uy) diverge with system size: 



dqxdqy q 




Gxxi<D 



■>pu{2TTaxLx/Ll) 



(4.14) 



if > 2TTaxLx; 
if Ll < 2TTaxLx 



where ^u{z) is a crossover function. The expression for 
(ity) is obtained by interchanging x and y in the equation 
for (u^). The anisotropy lengths ax and ay defined as: 



ax = [Kx/K)-- 



(4.15a) 
(4.15b) 



They diverge as the stability limit XxXy = X^y is ap- 
proached. 

The connected phonon correlation functions at two 
spatially separated points is 

CO (r) = ((^..(r) - u.-m){u,{v) - ^i,(0)))o (4.16) 
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For an infinite membrane 



-y'/{4a^\x\) 




2X'^ax y2y/ax\x\ J 

ifM«2y^, 
ifM»2v^. 



(4.17) 



wfiere erf(x) is the error function. The strong power- 
law growth of G% (r) indicates that thermal fluctuations 
lead to arbitrarily large relative displacements of two 
points distant on the membrane, again contrasting with 
the usual logarithmic growth with r in two-dimensional 
ordered xy-like systems. 

The stability of the spontaneous nematic (uniaxial) 
order can be analyzed by examining the rms fluctua- 
tions of the nematic director field ^n(r). Since in the 
nematic state, the director is "massively" tied to the 
antisymmetric part of the displacement gradient tensor 
rjij = djUi'l()\, orientational fluctuations can be com- 
puted from those of phonon fluctuations. To linear order 
in rotations of the director n away from its preferred ori- 
entation along the a:-axis, Suy — 9 — ^{dxUy — dyU^)- 
The angle correlation function is then 

Geeici) = \ WyGlM + ^'^^(q) - 2fe9,G|J,(q)] ■ 

(4.18) 

It is clear from Eqs. (|4.12a|l - H4.12cl) that within the har- 
monic approximation in-plane orientational fluctuations 
are finite. 

Turning to membrane out-of-plane fluctuations, we 
find that at harmonic level local rms undulations (/i(r)^)o 
behave the same as those of polymerized membranes with 

(Mr)^)o = 4 / / %^G",(q) 



27r J 22L 27r 



= 'J^^M^^l (4.19) 

where for simplicity we specialized the case of isotropic 
bending rigidity k and defined scaling function 



1 



with crossover property 



2n z^O 

2tt ^ , 



(4.20) 



(4.21) 



As in crystalline membranes, the strong growth 
implies instability of the flat phase to thermal fluctua- 



tions, as well as the importance of anharmonic elastic- 
ities, which (as in polymerized membranes) can stabi- 
lize the flat phase. Furthermore, in contrast to polymer- 
ized membranes, here the power-law divergent smectic- 
like in-plane phonon correlations that we found above, 
Eq. (|4.17|) . suggest that phonon elastic nonlinearities in 
Eq. H2.13|l may be important as well. We turn to these 
questions in the next two sections. 



V. ANHARMONICITIES AND THE 
e-EXPANSION 



As discussed in Sec. ^ rotational invariance in the 
target space requires the elastic free energy to be ex- 
pressed in terms of nonlinear rather than linear strains. 
The result is that there are anharmonic couplings in the 
elastic energy both among phonon fields and between 
phonon and height fields. The violent power-law fluc- 
tuations of these fields, the height field in particular, 
leads as in other membranes systems'^-T^ to diver- 
gences with system size L of perturbations in anharmonic 
couplings for membranes with spatial dimension D less 
than a critical value Dc- In this section, we consider 
the interplay between fluctuations and anharmonic cou- 
plings in _D-dimensional nematic elastomer membranes 
embedded in a d-dimensional space and show that per- 
turbation theory breaks down below Dc = 4. We then 
study the anomalous elasticity of these membranes in an 
e-expansion about D — 4:. Our interest is in develop- 
ing insight into properties of real two-dimensional mem- 
branes. We, therefore, consider only models that have 
a straightforward analytic continuation to D = 2. The 
two models we consider are the fully anisotropic Z?-axial 
elastomer, which has no surviving shear modulus, and 
the uniaxial model in which we set the shear modulus /i_L 
Eq. (|2.28() for shears in the plane perpendicular to the ne- 
matic direction to zero. These two models are equivalent 
in the physical limit of I? = 2. 



A. Perturbative analysis of elastic nonlinearities 

The elastic free energy Eq. (|2.37|) contains nonlin- 
earities associated with membrane undulations (involv- 
ing h field) and in-plane phonon nonlinearities. The im- 
portance of undulation nonlinearities is a consequence of 
membrane's vanishing tension (softness of out-of-plane 
undulations, controlled by curvature, rather than sur- 
face tension e nerg y) . As in an isotropic polymerized 
membrane [sL [T2| . undulation nonlinearities become rel- 
evant when the dimension of the reference space is lower 
than 4. To illustrate this point, we calculate the pertur- 
bative corrections to elastic constants Ay from undula- 
tion nonlinearities which are represented by diagram (a) 
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shown in Fig. ^ 

D 



D 



4-_D 



dcL 

2(4 -D) 



k,l 



d^D qlqf 



(5.1) 



where we assumed that the system has an equihbrium 
configuration of D dimensional sphere (disk in two di- 
mensional case) with radius L. dVLo is the differential 
surface element of the D dimensional unit sphere. Quite 
clearly, for D < 4, the fluctuation corrections 5\ij diverge 
with system size L. The associated nonlinear length 
scale, beyond which these corrections become compara- 
ble to Xij , and perturbation method breaks down is given 

by 



X 



(5.2) 



where A is the typical value of Xij and is defined as 
the angular integral in Eq. H5.1|) . 

If the modulus /i_L for shears in the direction perpen- 
dicular to the anisotropy axis in the uniaxial model is 
zero, then fluctuations in Ut Eq. 1)3. 6|l have the same q^* 
harmonic-theory divergence as height fluctuations, but 
with K{q) replacing K{q). Thus, when fi± = 0, both h 
and Ut contribute to divergences in Ay below D ^ A. 
This point was missed in the analysis of the fixed con- 
nectivity fluid fixed point in Ref. 6. 

As discussed in the Introduction, the spontaneously 
broken in-plane rotational symmetry of the nematic elas- 
tomer membrane leads to soft in-plane elasticity. As a 
result, in strong contrast to isotropic or crystalline mem- 
branes, in-plane nonlinearities in J^, Eq. H2.37|l . also cor- 
rect elastic moduli A,- 



Its contribution SXfj, as repre- 



sented by diagram (c) in Fig. ^ is given by 



D 



k,l 



(2^) 



\GlM)\ 



(5.3) 

In the Z?-axial model, the dominant fluctuations in Gj^^ 
are smectic-like, and it is straightforward to show that 
they cause the above correction to A"^ to diverge as 
Iji3-D)/2 i^qIq^ D = 3 if the phonon fluctuations retained 
their harmonic character down to D = 3. In contrast 
in the uniaxial (analytical continuation) model, the con- 
tributions to (5Aj" due to Uj fluctuations diverge below 
D = A, when ii± as discussed above. The nn part of 
C^j^ is, however, smectic-like and diverges more weakly 
as L(3-^)/2. 

Similar calculation shows that the perturbative correc- 
tion to in-plane elastic constants Kij is dominated by in- 
plane nonlinearities and also diverges below three dimen- 
sions. Thus the upper critical dimension for undulation 
nonlinearities is 4, while in-plane nonlinearities become 



relevant below D = 3. For 3 < D < 4, in-plane non- 
linearities are irrelevant, consequently in-plane curvature 
energy moduli Kij are only renormalized finitely. 

Undulation nonlinearities renormalize bending rigidi- 
ties Kij as well as in-plane elastic constants Xij. One 
would expect that as in the case of isotropic polymerized 
membranes, the perturbative corrections to also di- 
verge as L'^~^ . This is, however, not true, as one can 
see from the following argument. Note that if we neglect 
in-plane nonlinear terms (which is legitimate above 3 di- 
mension) in Eq. H2.37|l . and set Kij to be zero, we can 
integrate out the phonon fields completely and what is 
left is only the bending energy Hij{di^h ■ dj'^h) with- 
out any nonlinearities. Thus there should be no anoma- 
lous elasticity for bending rigidity k^ , were there no in- 
plane curvature rigidities Kij. This implies that, with the 
presence of ify , the renormalization of bending rigidi- 
ties Kij is dominated by anisotropic, smectic-like modes, 
for which in-plane curvature rigidities Kij are important, 
and the critical dimension below which Kij is infinitely 
renormalized is 3, as in smectic liquid crystals. 



(a) 




(b) 




(c) 



FIG. 1: Feynman diagrams renormalizing elastic constants 
A's (A), k's (B) and K^s (C). Solid lines denote undulation 
fields ft, dashed lines denote phonon fields u. 



The upshot of above discussion is that near D = 4 
in the D-axial model all in-plane nonlinearities in u are 
irrelevant, and in the uniaxial model nonlinearities in u„ 
and ul are irrelevant. Consequently, to capture the long- 
wavelength behavior we can use simplified expressions for 



11 



the full non-linear strain tensors given by 

D — axial : 

"ij ^{diUj + djUi + dih ■ djh), 

uniaxial : 
1 



(5.4a) 



{diUj + djUi + dih ■ djh + diUt ■ 9jUtjj5.4b) 

The effective Hamiltonian for studying membranes with- 
out a shear modulus is the Eq. 12.37|l with = and 
one of the above reduced strains. 

The difference between this model free energy and 
that of isotropic polymerized membraneslsHi^, as well as 
that of fixed-connectivity fluid membranes jSj, should be 
stressed. First, the bare elastic constants, Kij and \ij are 
anisotropic rather than isotropic, and the anisotropy in 
Xij cannot be eliminated by a simple rescaling of lengths. 
Second and more importantly, the energy cost for shear 
in the planes of anisotropy is zero because of the sponta- 
neous broken symmetry of the nematic state. Thus, our 
model free energy, Eq. (|2.37(l . is not a simple anisotrop- 
ically scaled generalization of fixed-connectivity fluid 6] . 
The matrix of coupling constants Ay is such the elastic 
energy can not be reduced to that of density variation 
alone. 



B. Renormalization Group and 4 — D-expansion 

In isotropic membranes below D = 4, height fluctu- 
ations lead to anomalous elasticity 0-^2 with bending 
modulus and elastic modulus, respectively, diverging and 
vanishing with wavenumber as 



K(q) 



A(q) ^ q 



(5.5) 



where the exponents rj^ and rjx are related via the Ward 
identity jg;] 



2?/k + ?7A 



D. 



(5.6) 



We have just argued that k is not renormalized for di- 
mensions near 4 in the models without shear moduli that 
we are considering. Thus 77^ = 0, and r]x ~ e. In this 
section, we will show explicitly within an RG calculation 
that this is indeed the case in both models we consider. 

We use standard momentum-shell renormalization- 
group procedures, integrating out a shell in momentum 
space with A/e' < g < A, where A is the ultraviolet cut- 
off, to produce fields u*^ (r) and (r) with wave-numbers 
with magnitude q < e~'A. We then rescale lengths and 
fields according to 



r'e'. 



h<{r) = e^'K'(r'), 



(5.7) 
(5.8) 
(5.9) 



so as to restore the ultraviolet cutoff to its original value 
A. In the uniaxial case, we decompose u as in Eq. (|3.t)|) 



and choose Xn = XL = X and Xt different from x- In 
the D-axial case, we can chose all of the Xi to be equal 
to X. It is convenient to choose the rescaling so as to 
preserve the nonlinear form of the strain Uij Eq. (|5.4|l . 
This requires 



X = 20 - 1 = 2xt 



1 



(5.10) 



The integration over the high wave- vector components 
of u and h can be performed perturbatively in nonlinear- 
ities of Hamiltonian, in a procedure very similar to that 
of perturbative analysis. The Feynman graph giving the 
1-loop corrections to Xij is shown in Fig. ^ diagram (a). 

After performing the rescaling and calculating the 
graphic corrections, we obtain the following RG flow 
equations: 

^ = (i?-2 + 2x)Ay - i^A,a4;Aij, (5.11a) 





-4H 




(D 


-4H 


^2x)Kp 


{D 


-4H 


K ^Xt)Klj 



(5.11b) 
(5.11c) 
(5. lid) 



fc,; 

nr 

dl 
dl 

where, for simplicity, we have set the ultra-violet cutoff 
A = 1. The components of the in-plane bending coeffi- 
cient, Kfp coupling to Un and ul scale with x, whereas 
those, , coupling to Ut scale with xt • The matrix M^i 
has similar but different forms for the uniaxial and D- 
axial models. In the D-axial model, Mki = dcM^i, and 
in the uniaxial model Mu = dcM^i + {D - 2)M^, where 



dn. 



= 



(27r)^ K^q) 
dflo qlqf 



(5.12a) 
(5.12b) 



{2Tr)D K^{q)- 

Note that Mj,;, like Xij and Aij defined in Sec. m is 
really a subset of the components constructed from the 
fourth-rank tensor 

dilo qtqjClkqi 



(5.13) 



{2^)D K^{q)qi 

and M^^i defined in a similar way. M^i does not trans- 
form like nor have the symmetries of a second-rank ten- 
sor. The Z)-axial and the uniaxial model differ mostly 
in their respective forms of Xij and . These different 
forms require slightly different fixed-point analysis, and 
we will now treat the two cases separately. 

Regardless of the model, we can choose (f) — {A — D)/2 
and Xt = (4 — D) /2 to keep and Kj.^ fixed. Then, 
X = 3 — D, and the inverse correlation function for the 
parts of u not in the anisotropy planes (i.e., not Uj) scales 
as 

G,r.i(q,A.„if,,) = e(^-s)'G,^-\e'q,A.,(Z),if,,(Z)) 

= e~''X,j{l)q,q,+K{q)q\ (5.14) 



12 



Thus, if Xij has a non-zero fixed point value A*,-, then, 
choosing e^q = 1, we have 



(5.15) 



Both Kij and Kl^ remain constant. 



C. The Uniaxial Model 

In the uniaxial model, both Eq. (|2.31() and My are 
uniaxial. Mij is easily calculated by taking the appropri- 
ate components of the full fourth-rank tensor Mijki Eq. 

My = Mii6aSji (5.16) 
+ (^^A/||^[(5,i(l - S,i) + (1 - Sn)Sji] 

The components of M^^ and M^^ of My have simple ex- 
pressions in terms of integrals over angle: 

dftjj cos"* 9 
dVtu cos^ 9 sir? 9 



^11 



Ml 



(Kin sin 9 



(5.17a) 
(5.17b) 
(5.17c) 



where is the angle q makes with the uniaxial axis (1- 
axis). A similar set of expressions applies to the compo- 
nents of with K{q) replacing K{q). Note that M23 = 
M^/[{D -l)iD + 1)] and M12 = M13 = M^^JiD - 1), 
whereas the components of the true second rank uniaxial 
tensor T,j satisfies T23 = T12 = T13 = 0. 

With these definitions, the recursion relations for the 
components of Ay become 



dXi 

~dr 

dl 

dh 
dl 



eAi - i(A?M|| + 2AiA2M||^ + A^Afj(J,.18a) 



eAa - -(A1A2M11 +A2M||^ 

+ AiA3A/||_l + A2A3M_l), 



(5.18b) 



eA3 - 2(^3^^^ + 2A2A3M||^ + XlM^ffi.lSc) 



where Ai, A2, and A3 were defined in Eq. H2.31|l . These 
equations have an unusual fixed-point structure as shown 
in Fig. 13 In appropriate two-dimensional planes in the 
three-dimensional space of Ai, A2, and A3, they exhibit 
the familiar four-fixed-point structure of systems with 
two coupled potentials in which there is one unstable 
Gaussian fixed point, one globally stable fixed point, and 
two fixed points that are stable in one direction and un- 
stable in the other. The full three-dimensional structure 



is topologically equivalent to what would be obtained if 
the two-dimensional structure is rotated about the axis 
connecting the Gaussian and globally stable fixed point 
and subsequently stretched anisotropically. In this pro- 
cess, the two mixed-stability fixed points become an el- 
liptical ring with stability exponent of zero for displace- 
ments along the ring and one positive and one negative 
exponent for displacements perpendicular to the ring. 
The fixed points and their stability exponents lui, 102, 
and IJJ3 are 



1. Gaussian 



Ai — A2 — A3 — 0, 
uji = UJ2 = ljJ?, = e. 



(5.19) 



2. Uniaxial nematic elastomer 



Ai - 



2eM_L 



A2 = — 



2eMi 



ll-L 



UJ2 



where /S.m 
3. Fixed Ring 

Ai = 



M||Mi 



Am 

CJ3 : 



A3 = 



2eM| 



II 



(5.20) 



2e 



M|| -|-2aM||_L +a'^Mx_' 
A2 = ckAi, A3 — a^Ai 
oJi — -e, W2 = 0, ^3 = 



(5.21) 



for —00 < a < 00. For every point on the ring, 
A1A3 = A2. T 
esting points: 



A1A3 = \\. The fixed ring includes various inter- 



(a) a — 1, Ai = A2 = A3. This is the fully 
isotropic fixed-connectivity-fluid fixed point of 
isotropic membranes 0. 

(b) a = -1, Ai = -A2 = A3. 

(c) a = ±00, Ai = A2 = 0, A3 = 2e/Mj_. 

(d) a = 0, Ai = 2e/M||, A2 = A3 = 0. 



D. The L»- Axial Model 

The analysis of the Z3-axial model is complicated by 
the fact that Ay has a large number {D{D + 2)/2 = 10 in 
13 = 4) of independent components. To simplify our pre- 
sentation, we will first consider flows in a restricted sub- 
space in which Ay is parameterized by only two parame- 
ters. Thus we will first find 4 fixed points with particular 
structure and show that one of them is globally stable 
and describes D-axial nematic elastomer membrane. We 
then give a general solution for all fixed points and show 
that it actually contains both the fixed point Eq. (|5.20l) 
and the fixed ring we found for uniaxial model. Further- 
more, we show that if anisotropy of Ay is turned on in 
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1. Gaussian 



FIG. 2: Renormalization group flows for the uniaxial model, 
showing the Gaussian (G) and uniaxial nematic elastomer 
(NE) flxed points and the fixed ring R, unstable in one direc- 
tion and marginally stable around its perimeter. The locus 
of flow lines, of which GA and GB are two examples, from G 
to R and extended beyond R form a distorted cone C. All 
points within and on the boundary of C flow to the stable 
fixed point NE. All other points flow to large coupling. 



the {D — l)-dimensional plane, the uniaxial fixed point 
Eq. H5.2(J|I become unstable and the system flows to the 
globally stable D-axial one. 

Our recursion relations are still given by Eq. (|5.11d|) 
without Kij. Choosing </> to keep Kij constant, we obtain 



dX 1 

^ = eA - -AAf A. 

dl = 2===' 



(5.22) 



where A and M are, respectively, matrices with entries 
Xij and Mij . As in the uniaxial case, this equation should 
have an isotropic fixed point will all Xij equal, i.e., with 
A ~ 1_ where 1 is the matrix with all entries equal to 
one. It is also clear that this equation has fixed point 
with A ^ M_~^ provided the symmetries of A and M_ are 
compatible. In the uniaxial case, M has more unequal 
components than A, thus the solution A ^ M_~^ is not 
permitted there. In the D-axial case, M has the same 
number of independent components as a D-axial second- 
rank tensor, and A ^ is permitted. We thus begin 
by seeking fixed points in the 2D subspace defined by 



A = XaM- 



Xhl. 



(5.23) 



The recursion relations for A^ and A^ are 
dX„. . 1 



- eAa - -Xl 



(5.24a) 

dX 1 ^ 

= eXb--Xh{2Xa + Xt J2 ^^^j)' (5-24b) 



where we used 1M_1 = "^-j Mijl. The fixed points and 
their stability exponents for these equations are 



Xa 
UJn 



Xb 

UJb 



2. Fixed-connectivity fluid 
Xa — 0, At 



2e 



UJb 



(5.25) 
(5.26) 



(5.27) 
(5.28) 



This fixed point is in fact the isotropic fixed- 
connectivity fluid flxed point found in Ref. (6| for 
isotropic membranes. 



3. Mixed 



Xa 



2e, 



2e 



oJa = -e, 



LUb = e. 



(5.29) 
(5.30) 



We call this flxed point mixed because its coupling 
constant matrix A has components characteristic 
of both the flxed connectivity fluid and the D-axial 
membrane fixed points. 



4. D-axial nematic elastomer 

Xa = 2e, Afc = 0, 
LUa = u;b = -e- 



(5.31) 
(5.32) 



This fixed point is in fact globally stable, i.e., it 
is stable in all directions. Plugging A — XnM~^ 
into the recursion relation Eq. (|5.22|) and lineariz- 
ing, we obtain dSX/dl = —eSX. Thus, the stability 
exponent is — e for all directions. 

We have thus identified one globally stable fixed point, 
the nematic elastomer fixed point, and two others. The 
full fixed point structure of the flow equation Eq. 1)5.221) 
is actually not difficult to determine. Since M is generi- 
cally positive definite, we can define a new "(symmetric) 
coupling constant matrix P — MyXMy / (2e), whose flow 
equation given by: 



dl-- 



-P = P - 



(5.33) 



One immediately see that every projection matrix(3J| is 
a flxed point for P_ and vice versa. This means that the 
general solution for the flow equation Eq. 15.22)1 is given 
by: 



A = 2eM-5PM~5. 



(5.34) 



with P an arbitrary projection matrix. 

We can classify fixed points, or more generally fixed 
subspaces, by the dimension Dp of the space that P 
projects onto. Of course, TrP — Dp. If the dimension D 
of the elastic manifold is an integer, then P can project 
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onto all subspaces with dimensions Dp = 0,1, 
When D is not an integer, the classification is less clear. 
A convenient set, however, is the set with dimensionali- 
ties Dp = 0,1, [D], and Dp = D,D - 1, ...,D- [D], 
where [D] is the greatest integer less than or equal to 
D. li Dp = D, = / is the unit matrix, where the 
superscript indicates the dimension of P. If Dp — 1, 
= CiCj for any unit vector e. A (Z? — l)-dimcnsional 



operator P-j^ 



can also be constructed from 



the unit vector e. Similarly, k- and {D — fc)-dimensional 
projection operators can be defined via Pjj ~ J2i=i ^u^ij 

and Pj;f~''^ = Sij - P^j, where ■ e; = Ski- 
To study stability of the flow equation, Eq. (|5.33|) , for 
P for a given fixed-point projection matrix, P^° , we ex- 
press deviations of P from PP'-' as 



6P = P-Pt 



P. +P„ 



(5.35) 



the space of coupling constants A. This space is neces- 
sarily compact since the subspaces defined by P^" are 
parameterized by unit vectors. 

There is only one projection operator with Dq = D. 
This is the operator P^ = / that projects onto the whole 
space. For this case, P2 and P3 are both zero, and the 

fixed point is stable in all directions. It is the globally 
stable fixed point with A = 2eM~^, which is identical to 
the stable D-axial fixed point of the restricted set class of 
couplings defined by Eq. (|5.23() . The other fixed points 
for the restricted set of couplings must correspond to 
some P^° with Dq < D. It is straightforward to show 
that the fixed-connectivity-fluid fixed point corresponds 



with 



El 



(5.42) 



where 



^1 



^2 

P„ 



=0 "==0 ' 

a- E^°)sm- E^°) 



pDo 



)SPP. 



Do 



(5.36) 
(5.37) 
(5.38) 



Recall that P is a _D-dimensional symmetric matrix with 
D{D + l)/2 independent components. P^ is the pro- 
jection of 6P onto the Z3o-dimensional subspace defined 
by and it has Do{Dq + l)/2 independent compo- 

nents. Similarly, P^ is the projection of SP onto the 



D 



Dq dimensional subspace defined by / — P^ " with 



(D — Dq){D — Dq + 1)/2 independent components. Fi- 
nally, P^ represents the Dq{D — Dq) components of 5P_ 



that couple the subspaces defined by P^° and / — P 



To linear order, the flow equations 



for (5°P 



are 



±P 
Ap 

AP 
d/=3 



^2' 



= 0. 



(5.39) 
(5.40) 
(5.41) 



Thus, P^ is stable, P^ is unstable, and P^ is marginally 
stable. This means that the I^o-dimensional fixed-point 
projection matrix is stable with respect to any change 
in P in the Dq{Dq + l)/2 dimensional space of Dq x Dq 

matrices that operate in the space defined by P^"; it 
is unstable with respect to any changes in P in the 
(£'--Do)(-D--Do + l)/2-dimensional space of {D - Dq) x 
{D — Dq) matrices that operate in defined by / — P^"; 
and it is marginally stable with respect to changes in P 
in the Dq{D — Z?o)-dimensional space of matrices that 
couple P^° to 6 — P^° . These results imply that the set 
of fixed points defined by all Z3o"dimensional projections 
matrices is a Z3o(Z?— Z?o)-dimensional surface in the space 
of all possible symmetric matrices P or, equivalently in 



where M^^- is the ij component of the matrix . Thus, 
this fixed point is actually a single point in a (13 — 1) = 
(3 — e)-dimensional fixed manifold. Similarly, the mixed 



fixed point corresponds to Pr 



D-l 



j — CiCj. There are 
other unstable fixed points for A not described by the 
restricted set defined by Eq. H5.23|l . in particular those 
with Do = 2 or £) - 2. 

We have identified all of the fixed point manifolds of 
the D-dimensional coupling matrix A. These must in- 
clude the fixed points of the uniaxial model discussed in 
the preceding subsection. When uniaxial constraints are 
applied, it is natural to construct unit vectors and pro- 
jection matrices from the components of M parallel and 
perpendicular to the anisotropy axis. It is not difficult to 
show that the stable uniaxial fixed point corresponds to 
a two-dimensional projection operator 



where 



-Poij ^ eiiEij + e2ie2j, 



Mu 1* 



(5.43) 

(5.44) 
(5.45) 



Thus, the uniaxial fixed point is a point, satisfying uni- 
axial constraints, in a 2{D — 2)-dimensional fixed man- 
ifold of all possible couplings. A point on the uniax- 
ial fixed ring parameterized by a corresponds to a one- 
dimensional projection operator P^ defined by the unit 
vector 



lia) 



ik 



(M|| -t- 2aM||j^ -f a2Af_L) 



(5.46) 



The set of vectors e{a) defined by all a define a one- 
dimensional loop in a (13 — l)-dimensional fixed manifold 
in the space of all possible D-dimensional couplings. 
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VI. DISCUSSION AND CONCLUSION 



I T 




FIG. 3: A possible phase diagram for ideal nematic elastomer 
membranes. As the temperature is lowered a crumpled mem- 
brane undergoes a transition to isotropic flat phase at Tcf, 
followed by a 2D in-plane isotropic-nematic like transition 
to an anisotropic(nematic) flat phase. As T is lowered fur- 
ther, this anisotropic flat phase becomes unstable to a ne- 
matic tubule phase, where it continuously crumples in one 
direction but remains extended in the other. At even lower 
temperature, a tubule-flat transition takes place at Ttf- 



In this paper we studied thermal fluctuation and non- 
linear elasticity of nematically-ordered elastomer mem- 
branes in their flat phase. For the physical case of 
two-dimensional membranes, we found that at harmonic 
level in-plane phonon correlations are short-ranged but 
the in-plane orientational order remains long-range, in 
spite of violent thermal fluctuations. A generalization 
of a nematic elastomer membrane to an arbitrary D- 
dimensional membrane with either uniaxial or Z?-axial 
nematic order, allowed us to study the effects of out-of- 
plane undulation and in-plane phonon nonlinearities. We 
found that undulation nonlinearities are relevant D < 4 
and dominate over the in-plane nonlinearities that only 



become important D < 3. Focusing on the dominant 
undulation nonlinearities (and neglecting in-plane non- 
linearities) we performed an RG calculation combined 
with an expansion about _D = 4 and found that for 
3 < Z? < 4, bending rigidities are only finitely renormal- 
ized, while in-plane elastic moduli become singular func- 
tions of a wavevector (i.e., exhibit anomalous elasticity), 
vanishing with a universal power-law. This power-law is 
controlled by an infrared stable fixed point whose sta- 
bility we analyzed in detail for the uniaxial and _D-axial 
analytic continuations, finding agreement in fixed point 
structure. This analysis also allowed us to make contact 
and recover some of the results previously obtained in 
the studies of (crystalline) polymerized membranes. In 
particular we found that the so-called connected fiuid is 
realized as a fixed point of a nematically-ordered elas- 
tomer membrane that is unstable to the globally stable 
nematic-elastomer fixed point. 

Despite of some of the success in our understanding of 
the behavior of nematic elastomer membranes, there are 
obvious limitations of our analysis, most notably in the 
application of our work to the physical case oi D = 2 elas- 
tomer membranes. This shortcoming primarily has to do 
with the neglect of in-plane elastic nonlinearities, which, 
near the Gaussian fixed point become relevant for D < 3. 
While it is very likely that the subdominance of these in- 
plane nonlinearities relative to the undulation ones will 
persist some amount below D = 3 |32j | , we expect that in 
the physical case of D = 2 all three nonlinearities need to 
be treated on equal footing. Doing this remains an open 
and challenging problem. 

We conclude with a discussion of the global conforma- 
tional phase behavior of nematic elastomer membranes. 
As with ordinary polymerized membranes we expect, 
upon cooling, isotropic elastomer membranes to undergo 
a crumpling (flattening) transition from the crumpled 
to flat-isotropic phase. Upon further cooling, an in- 
plane (flat) isotropic to (flat) nematic transition can 
take place. As shown by Toner and one of us 23], poly- 
merized membranes with arbitrary small amount of in- 
plane anisotropy inevitably exhibit the so-called tubule 
phase whose properties and location in the phase diagram 
are intermediate between the high-temperature crumpled 
and low-temperature flat phases. Thus we expect that 
for nematically-ordered elastomer membranes there is a 
similar nematically-ordered tubule phase. Since in such 
a state the in-plane rotation symmetry is spontaneously 
(as opposed to explicitly) broken, we expect qualitatively 
distinct in-plane elasticity distinguished by the presence 
of a new in-plane soft phonon mode. Consequently a 
nematic tubule should be qualitatively distinct phase of 
elastic membranes. This discussion is summarized by a 
schematic phase diagram for a nematic elastomer mem- 
branes, illustrated in Fig|2| A detailed analysis of the ne- 
matic tubule phase and phase transitions between it and 
other phases will appear in a separate publication j25j| . 
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APPENDIX A: SCALING OF HARMONIC 
FLUCTUATIONS AT 2 DIMENSION 

In this appendix, we show detailed calculation to jus- 
tify the approximate form of harmonic phonon propaga- 
tors in 2 dimensions, as presented in Sec. 14.121 

The propagators G9 (q) are easily found through 
equipartition or by a Gaussian integration 



GL(q) = 



1 



,<I>o(^,^), 



1 



>3x 



y 1^- % 

^ ,o.xql ayql 
$o( -,^^). 



Xxij4>q{ 



-M 



axql ayql^ 



qx 



Qy 



<iv 



{\xql+Kxql){\yql + Kyqty 
with anisotropy lengths ax and ay defined as: 
(i^./A,)5, 

-■{KylXy)'^. 



ax 
ay 



(A2a) 
(A2b) 



Because of the anisotropy of the nematic state the prop- 
agators are highly nontrivial even at the harmonic level. 
Their angular dependence is encoded by the crossover 
functions 



00(2) = 



1 



1-l-z^ 



1 - P^0O (2^)00 (w) 



(A3a) 
(A3b) 



illustrated in Fig0]and with asymptotics of the double- 
crossover function ^q{z^ w) given by 



with the ratio 



1 



if z 
if z 



and w 
00 or w - 



> 
00. 



\2 

'^xy 

\x\ 



(A4) 



(A5) 





1 ^^^X 







qx 



FIG. 4: A polar plot of the crossover function 
^o{axql/qx,ayq^/qy) with q = i/ii4~gf fixed. Shaded 
are the "decoupled region" where \qx \ <^ a^qy or \qy \ <^ ayq^. 
While $0 is positive and finite for all q, in the limit 5 — > it 
exhibits cusps at On = SiTctan{qy / q^) = nm 12 (n = 0, 1, 2, 3). 



required by stability to be less than 1. 

As can be seen from the asymptotic form given in Eq. 
(|A4|) . the scaling function ^o{axql/qx,o,yq'^/qy) is finite 
for all q and therefore simply provides an angular mod- 
ulation to the phonon correlation functions G^j{q). Its 
value ranges between the "decoupled" and "coupled" val- 
ues of 1 and 1/(1-/9^), with decoupled regime defined 
by a union of \qy\ <C a^g^ and \qx\ ^ a^q^ regions in 

q. The coupled regime is the complement of the de- 

— - 1 

coupled regime Uxq^ ^ \qx\ ^ ay^\qy\^, as illustrated 
in Fig^ As a result, seZfcorrelation (diagonal) func- 
tions G^j. (q) and Gyy (q) are essentially those of two in- 
dependent 2D smectics with x— and y— directed layer- 
normals and corresponding phonons Ux and Uy, respec- 
tively. The only effect of the cross-coupling Xxy on these 
phonon seZ/-correlation functions is to finitely enhance 
their amplitude in the coupled regime, without modify- 
ing their long-wavelength pole structure. Thus in order 
to study the fluctuation of Ux field, we only have to con- 
centrate in wave vector region q^ ^ 0,'^qy- For <C a^fZy, 
^o{axql/qx,ayql/qy) « 1, and G°(q) « l/Kxqy, while 
for \qy\ -^ql-^ alq^, <S>o{axql/qx,ayql/qy) « l/(l-p2), 
and G2(q) ~ l/A'^q^. This is exactly the same as the 
asymptotic behaviors of Eq. (|4.12a() . Thus we see that 
Eq. (|4.12a|) is a good approximation for G'^ in wave vec- 
tor region q^ ^ Of course for q^ S> g^, the ratio 
between Eq. f l4.12a|) and GjJ is approximately 1/(1 — p^), 
but this region is not important for fluctuation of Ux 
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phonon anyway. This analysis also applies to G ° if we 
exchange label x with y in every place. Thus Eq. f l4.12bp 
is also a good approximation for G°. 

In contrast, the phonon cross-correlation (off-diagonal) 
function G^y (q) depends strongly on whether q is in de- 
coupled (union of <C CLyq^ and \qx\ "C axq^ regions) 

or coupled {a^Qy <C \qx\ ^ \qylay\^) regimes. At long 
scales, near the two dominant x and y smectic regions of 
q, it is strongly subdominant to the self-correlation func- 
tions G'°j.(q) and G°y(q), down by a factor of qy and g^., 
respectively. 

The subdominance of the cross-correlations relative to 
the self-correlations can also be seen by analyzing the 
behavior of the cross-correlation ratio of Ux and Uy 



^GO,(q)GO,(q) 
Simple analysis shows that in the decoupled regime 



(A6) 



Pxy(q) 



«P, for \qx\ ^axql 



(A7) 



suggesting that thermal fluctuations of Ux and Uy are 
nearly independent. On the other hand, in the coupled 
regime 



Pxyil) ^ P, for ttxql < < ky/ayl- 



(A8) 



and Uy are strongly correlated. In this region, we have 



^0 ^ 



(A9) 
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